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ABSTRACT 


One transitional step for the development of a 1 MW power directed energy 
weapon is the proposed 100 kW upgrade of the Thomas Jefferson National Accelerator 
Facility’s Free Electron Laser (FEE). To improve the performance of the EEL, the use of 
the step-taper undulator is explored. Steady-state gain, final steady-state power, and the 
induced electron spread as a function of desynchronism and taper rates are determined. 
Comparisons are made to the conventional periodic and linearly tapered undulators. The 
multimode simulations used showed that the TJNAE lOOkW EEL is feasible. Simulations 
results with Q =10 show that the inverse step-taper undulator A = -7i achieved the 
highest final power of 190 kW at a desynchronism value of d = 0.01, while maintaining 
the induced energy spread well below the engineering limit. The validity of our results is 
verified against experiments conducted in the TJNAE EEL facility. The simulations and 
the experimental data are in good agreement and consistent with analytic theory. 
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EXECUTIVE SUMMARY 


The Free Electron Laser (EEL) weapon is excellent for providing state-of-the-art 
air defense. It is also effective, if need be, in offensive operations. It is possible for this 
weapon system to achieve a soft kill whenever needed by adjusting lethality. When 
combined with its accuracy makes it a valuable tool in operations where the requirement 
is not to only win, but also to avoid casualties. 

One transitional step toward the 1 MW output power required for a laser weapon 
is the proposed 100 kW upgrade of the Thomas Jefferson National Accelerator Eacility’s 
EEL. At present, the TJNAE EEL is the most powerful EEL in the world and operates at 
more than 2 kW average power. 

Our multimode simulations showed that this development is feasible, and for the 
current TJNAE lOOkW EEL parameters with 2 = 4.2, the inversely tapered undulator 
A = -71 achieved the highest final power at a desynchronism of d = 0.04. The output 
power produced was P -120 kW, which is above the lOOkW objective. The small values 
of the induced electron energy spread (-3.5%) made it possible to explore a design with 
lower 10% losses of (2 =10) in an optical cavity. The higher Q factor was chosen to 
increase the efficiency of the system, while maintaining the induced energy spread well 
below the 15% engineering upper limit. Simulations results with 2 =10 design indicate 
that the inverse step-taper undulator A = -71 achieved the highest final power 
P = 190 kW at a desynchronism of d = 0.01. The inverse linear taper 6 = -271 achieved 
the highest final power P =200kW at a desynchronism of d=0.01. Eor all the 
aforementioned cases the steady-state gain, the final steady-state power, and the induced 
electron spread as a function of desynchronism and taper rates were determined. 

The validity of our results was verified against experiments conducted in the 
TJNAE EEL facility at lower power. EEL operation as a function of various taper rates 
and desynchronism values was studied. The purpose of the simulations was to compare 
experimental observations with theoretical analysis, and, by using a wider range of 

xvii 



parameters than allowed in the experiment, to extract more physical meaning from the 
results. 


The simulations and the experimental data agreed well and for the most part were 
consistent with analytic theory. This theory indicates that as the taper increases, the 
desynchronism width decreases. The experiment shows this trend, but not as clearly as 
the simulation data. Only one point shows a significant deviation between the simulation 
and experiment. This has resulted in more examination of the experimental data, 
revealing that the conditions changed during the measurement of the point in question. 

The maximum gain with no tapering, 5=0, occurred, as would be expected. As 

the taper rate 5 increases or decreases, the gain decreases from 160% to about 40% 
symmetrically around 5=0. At positive or negative taper rates near 5 = ±871 , the gain 
plateaus observed decreased much more slowly with an increasing taper rate. As the taper 
rate increases in magnitude from |51 = 871 out to |51 = 2471 , the gain decreases from about 

40% down to about 10%. The gain plateau at large values of taper begins at a value of 
taper 5 that causes the gain spectrum to change shape so that there are two peaks with 
nearly the same peak gain. As the taper is increased and the peak available gain 
decreases, the gain spectrum acquires multiple peaks with comparable gain. As the 
primary gain peak decreases, other surrounding gain peaks increase in comparison. This 
can also occur in the untapered case, but results in a smaller gain peak. As the peak 
available gain increases, or decreases, at any value of tapering, the desynchronism curve 
width Ad increases or decreases correspondingly. With more gain available there are 
more values of desynchronism that are above the threshold which makes the 
desynchronism curve wider. 
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1. INTRODUCTION 


Directed energy weapons (DEW) is not a new initiative that “just appeared.” It is 
the consequence of changes in political philosophy and technology. The emergence of 
more and more power sources of microwave energy is clearly a key factor in their 
applicability. Furthermore, many of these sources are available at reasonable cost and 
satisfactory sizes. 

DEWs have been on the on the U.S. military’s “love to have” list for several 
decades. Within the United States, all three services have been actively developing laser 
technology. One of the most impressive land-based Tactical High-Energy Easers (THEE) 
is the “Nautilus”. It is a joint U.S. Army and Israeli project, which uses a Mid-Infra-Red 
Advanced Chemical Easer (MIRAGE). At present, it is being tested at White Sands 
Missile Range in New Mexico and has been successfully tested against the supersonic 
BM-21 122mm artillery rocket. The Nautilus program is still being tested and reached a 
milestone when, in June 2000, a test firing resulted in the destruction of a Katyusha type 
missile in flight [Ref. 1]. 

The U.S. Air Force, as well as the U.S. Army, in conjunction with the Air Borne 
Easer (ABE) program, adapted the chemical-oxygen-iodine laser (COIE) to a flyable 
model against missiles. This project, which was on schedule and within budget, is now 
transitioning from the early concept design to operational use. 

The U.S. Navy is also interested, and on the 24th of April 2001 the Commander in 
Chief of the U.S. Atlantic Fleet sent the following letter to the Chief of Naval Operations: 

...and has proven that Speed of Eight weapons can be very effective 
against these small high speed threats. Such a laser weapon would offer 
our Naval forces an extremely versatile weapon to counter numerous soft 
and hard targets. A High Energy Easer weapon can be designed to deliver 
energy that can track, warn, damage, mission kill, and if need be, destroy a 
threat. I believe it is exactly this type of weapon system that our forces 
need in the littoral environment where, even though the threat may not 
always be as sophisticated as a highly maneuverable cruise missile, 
intentions are often more difficult to discern and timelines are extremely 
short. 
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A Free Electron Laser (EEL) as a weapon is an excellent approach to fulfilling 
present and future requirements. It could provide state-of-the-art air defense not only to 
ships but also to sensitive areas such as Naval and Air Eorce Bases. It is capable of 
achieving a high-probability hard kill against all cruise missiles. It increases the 
engagement range and minimizes engagement time and thus reduces the threat of 
simultaneously arriving cruise missiles. Reliability is almost 100%. Since only light and 
vacuum are involved in the weapon system, minimal maintenance is required. Its infinite 
magazine, its instantaneous reaction at the speed of light time and the rapid re¬ 
engagement point offer many advantages over gun and missile-based close-in weapons. 

It is also possible for this weapon system to achieve a soft kill whenever needed 
by adjusting the lethality. The latter, combined with its pinpoint accuracy, makes it a 
valuable tool in operations where the requirement is not to only win, but also to avoid 
casualties [Ref. 2]. This fact is present in all peacekeeping operations, as it was in 
Kosovo and now in Afghanistan. This philosophy has been another push toward the 
deployment of DEWs and EELs. 

Presently, EEL research is widespread and advancing on many fronts. The 
impetus comes from some unique advantages the device offers over other types of lasers. 
[Ref. 3]: 

• “Potential for producing very high power by extending technology 
developed for existing electron accelerators. The evolution of power 
densities in magnetrons, klystrons and gyrotrons have reached their limits 
while the limit for EELs have not yet been reached as Eigure 1 indicates. 
There the product of the average power and the square of the radiation 
frequency, as a figure of merit, is shown. These available microwave 
powers were unthinkable just a few years ago. 

• Prospects for wall-plug efficiency more than 10 percent or more at high 
powers. 

• Potential for tunable operation from the millimeter-wave region to the 
extreme-ultraviolet or soft x-ray region, although no single device would 
operate over such a broad range. (The broad tuning range comes from the 
fact that the electrons are not bound when they emit light; it would be 
achieved by varying the electron energy and the spacing and strength of 
the magnetic field)" 
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Figure 1. Progress in the Development of a Variety of Vacuum Electronic Sources of 

Coherent Radiation as Measured by the Evolution of the Product of the Average Power 
and the Square of the Erequency. “Erom [Ref. 4].” 

Chapter II of this thesis describes the EEE weapon concept and the improvements 
gained in the field. Present limitations mentioned at the Workshop 2001 in Virginia, 
USA, compare the PEE with other weapons (HEE and conventional) and address its 
benefits. 

Chapter III gives an overview of the theoretical background and discusses the 
physics behind the PEE. 

Chapter IV describes the results of simulations for the proposed 100 kW upgrade 
to the TJNAE PEE. This upgrade is a transitional step for the 1 MW output power. This 
work was originally presented at the 23rd International PEE Conference in Darmstadt, 
Germany and has been accepted for publication in Nuclear Instruments and Methods in 
Physics Research. 

Chapter V presents the analysis of simulation results based on experiments 
conducted at TJNAE. Simulations and experimental results are compared to analytical 
theory. 
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n. FEL 


The FEL is a device which extracts light energy from a beam of relativistic 
electrons passing through a spatially periodic magnetic field. John Madey first introduced 
the concept in 1970 [Ref. 5]. It is unl ik e other lasers, such as solid-state, semiconductors 
and liquid lasers, because it relies cn light emitted by electrons that are not bound to 
atoms. 

An FEL consists of an electron beam source, a periodic transverse magnetic field 
(an “undulator” or “wiggler” magnet), and an optical resonator as shown in Figure 2. 

Electtx)n beam 



resonator 



mirror 


Figure 2. The FEL Schematic. 


The undulator magnet in an FEL imposes a transverse acceleration on the 
electrons, resulting in the generation and amplification of light, which co-propagates with 
the electron beam through the magnet. As in a conventional laser, the resonator mirrors 
provide feedback around the amplifying medium, permitting the stored radiation to build 
up to saturation and produce a very powerful, coherent laser beam. 

A. WEAPON CONCEPT-PRINCIPLES INVOLVED 
I. The Way the Target is Destroyed 

A laser beam in the basic fundamental mode closely approximates a Gaussian. 
The intensity drops off from a maximum and has no side lobes. The typical laser beam is 
quite narrow at no more than a few millimeters in diameter. Furthermore, when operating 
at that mode, the minimum waist of the optical beam is where the highest energy density 
occurs. 
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Since the beam resembles a truncated plane wave, it is, of course, highly spatially 
coherent. Its directionality may be thought of as a manifestation of that coherence. Laser 
light is quasimonochromatic, and generally has an exceedingly narrow frequency 
bandwidth indicating that it is highly temporally coherent. Another attribute is the 
amount of radiant power that can be delivered in that narrow frequency band. The laser is 
distinctive in that it emits all its energy in the form of a narrow beam. In contrast, a lOOW 
incandescent light bulb may give out considerably more radiant energy than a lower- 
power cw laser, but the emission is incoherent, spread over a large solid angle with a 
broad bandwidth. 

A good lens can totally intercept a laser beam and focus all of its energy into a 
minute spot whose diameter varies directly with wavelength X . Spot diameters of just a 
few thousandths of an inch can readily be attained with lenses that have a conveniently 
short focal length. The laser intensity that can readily be generated in a focused beam is 
over 10^^ W/crn^, in contrast to an oxyacetylene flame having an intensity of roughly 10^ 
W/cm . A focused laser beam of a few kilowatts can bum a hole through a quarter- inch 
stainless steel plate in about 10 seconds. [Ref. 6] 

Studies and scaling experiments have shown that a liter of a missile’s material can 
be destroyed in a few seconds with 1 MW of laser power [Ref. 7]. One liter corresponds 
to a hole with the dimension 10cm x 10cm x 10cm which is generally enough to 
structurally disable the missile. Missile destmction can be achieved in different by an 
explosion of the missile warhead, or an aerodynamic instability causing the missile to 
breakup. 

A targeting system on a launching platform will control the beam directed to the 
missile and bum a hole in it. Unlike conventional bullets or missiles, light instantly 
travels to the target and does not suffer from gravitational effects. The Acquisition 
Tracking and Pointing (ATP) system is in use and has been successfully tested in the 
MIRACL project many years ago [Ref. 1]. Acquisition includes the detection of the 
target in the tracking system using infrared radiation. Tracking begins when a series of 
consecutive observed positions collected to allow for the filing of the azimuth, elevation 
and time. Connecting these points permits the constmction of nominal track data that 
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point the telescope in order to keep the image of the target centered in the tracking focal 
plane. The system temporal response, which includes a camera processor and motors, 
must be fast enough so that the target does not disappear from the camera field of view. 
The pointer/tracker will be stabilized against unintended motion by using a high speed 
gyroscope and inertial measurement system. ATP could cost $6 million, including the 
multiple target acquisition sensor and stabilization system [Ref. 8]. 

2. The Threat - Requirements from a Weapon 

Targets are evaluated with respect to their maneuverability, size, and vulnerability 
to attack. Target maneuverability imposes a requirement that the missile must be capable 
of terminal maneuverability to effect lethality. 

The airborne-target spectrum is quite varied and consists of low altitude, slow- 
moving targets in ground clutter, high-altitude fast-moving targets, head-on targets with 
tremendous closing speeds and side-encounter targets with high-g turning requirements. 
If the threat is other than airborne-targets, no known limitations are implied. In contrast, 
the surgical accuracy of the system, the absence of fragments and the instant battle 
assessment make the system operable for many targets and environments. Range and 
time-to-range is a factor that does not greatly affect the PEL weapon design. The 
engineer need not consider design features affecting range and time-to-range such as 
drag, wings, controls, lifting surfaces, and other parasitic losses. The incredibly small 
time of flight of 33 |isec, or the time it takes the beam to reach the missile at 10 km, is the 
key issue. However, environmental constraints and atmospheric transmittance have to be 
taken into account. 

Cost per shot must be much less than cost of threat target, with high lethality per 
shot. For example, an Unmanned Aerial Vehicle (UAV), which is used for 
reconnaissance, may be worth less than an expensive sea-sparrow missile. The 
motivation for cheap, effective defense does not imply that we are not going to use all the 
means of defense, whenever human lives are placed in danger. Additionally, the 
individual engagement time must be short in order to engage a large number of targets. A 
dwell time of 3 seconds makes it possible to handle several incoming supersonic targets 
simultaneously. 
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3. FEL Comparison with Existing Weapons 

Present systems detection is limited to the horizon for low level flying threats. 
Unlike the antiballistic-missile kinetic kill systems designed for missile-terminal phase- 
defense, the FEL mode of attack offers several key advantages. First, the target is much 
slower moving and there is less urgency for targeteers to worry about distinguishing 
decoys. Additionally, a significant advantage of the FEL is that it is environmentally 
friendly as it does not generate the toxic effluents common to other types of weapons. 
This issue of “the exhaust gas contaminants” is serious since the crew could be harmed, 
and equipment damaged. Furthermore, another disadvantage is the need for a logistics 
supply of hazardous and toxic chemicals. The latter probably is the strongest 
disadvantage since no one today wants to handle such materials. 

With pinpoint accuracy, FELs can destroy precise targets in rapid order, while not 
harming non-military facilities and people. This is a benefit when compared with other 
weapons, which often destroy more than the intended target. Thus, the damage to 
surrounding populations and facilities is minimized. 

With an almost zero reaction time, FELs enhance the survivability of the user, 
especially in cases when the threat is unexpected. This reaction time at long range 
disables the enemy threat before getting near. Current close-in weapons, such as the 
Phalanx gun, suffer from dispersion as a result of vibrations. Moreover, simulations and 
demonstrations have shown that the probability of hitting the missile reduces as range 
increases. The Phalanx does not destroy the missile at a distance sufficiently far enough 
away to protect the ship and the crew from serious damage [Ref. 9]. 

The Rapier, used in the 1982 Falklands campaign, the Hawk, and the Arrow 
systems are deployed worldwide. They can provide short-range air defense of airfields 
and smaller troop concentrations. These systems are more than 20 years old, and even the 
famous Patriot air defense system is not as good as thought to be. “Patriot wasn't what I 
would call highly effective. We lost some lives to ballistic missiles during Desert Storm. 
Twenty-six Army soldiers were killed in the barracks in Dhahran as a result as an attack 
from a Scud missile” as stated by Ft. Gen. Lester L. Lyles [Ref. 10]. 


8 



FEL is designed in such a way that maintenance is minimal. It can be fired 
instantly and no other maintenance is needed. Furthermore, it could play a role as a 
“weapon of choice” in today’s U.S. military because it has the ability to set power levels 
on a graduate scale and provide the user with choices unavailable in kinetic energy 
weapons as well the ability to minimize casualties by permitting a selection of non-lethal 
options. It has speed-of-light engagement with a variety of targets and has the potential to 
produce a range of precisely controlled effects, as well as the potential for deep 
magazines and low cost per shot. The absence of costly rounds or defensive intercept 
missiles causes the cost per shot to drop to ~ $2/missile (using Igal of fuel), while the 
RAM is ~$0.7 millions/missile and the PHALANX is ~$3k/missile. 

B. FEL DESIGN 

1. Application 

The application requires high power in a small footprint. The proposed FEL 
consists of six basic components: the electron injector, the linear accelerator, the 
undulator, the resonator, the refrigerator and the electron beam dump. All these 
components can be packaged into a box with dimensions of 12m x 4m x 4m as shown in 
Figure 3. 


LASER BEAM 


4 m 



REFRIGERATOR 




12 m 



FEE MW Design. (After [Ref. 11]). 
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Figure 3. 




















The configuration and operational parameters are provided in Table 1. 


Electron beam energy 

E^=140 MeV 

Relativistic Lorentz factor 

Y=275 

Average current 

I = 0.6 Amp 

Accelerator RE frequency 

a =750 MHz 

Accelerator gradient 

17 MeV/m 

Peak current 

7 = 800 Amp 

Charge 

q = 0.8 nC 

Electron beam radius 

Ve = 0.08 mm 

Electron beam density/pulse 

p= 8xl0'Vcm^ 

Optical mode waist radius 

Wo = 0.08 mm 

Number of undulator periods 

N = 20 

Undulator period 

\ = 3 cm 

Undulator length 

L= 60 cm 

Undulator gap 

g = 1 cm 

Undulator parameter 

K = 2 

Resonator length 

5=12m 

Optical wavelength 

A- =1 pm 

Quality factor 

2 = 4 

Extraction efficiency 

^ ==1/4A = 1.25% 

Energy spread 

Ay /y ~ 2//V =0.1% 

Rayleigh length 

Zo = 1.8cm 

Intensity on mirrors 

2l0kW/cn? 

System dimensions 

12mx4mx4m 

Weight 

60 Tons 

Cost 

$60 million 


Table 1. MW PEL Parameters. (From [Ref. 11]) 
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Some of the challenging issues are intensity in the mirrors, effective recirculation 
of the 84 MW electron beam and beam transportation. A deployable MW-class PEL will 
need a sophisticated feedback, control and diagnostic system. Due to the short turn-on 
time and limited mission duration, the system must be self-diagnosing, and self- 
correcting for shock and vibrationally induced misalignments. [Ref. 12] 

2. Atmospheric Propagation 

Although powerful at the point of origin and delivered at more than 186,000 miles 
per second, a laser’s energy can be absorbed and diffused even in a relatively short 
passage through the atmosphere. There are propagation issues that need to be considered: 

• Thermal blooming, 

• Absorption, 

• Diffraction, 

• Scattering. 

It is possible to overcome these difficulties by the use of scaling laws and various 
simulations to predict what will happen in the MW class. These tools have limited 
applications as scaling and simulations do not always provide reliable results. The use of 
adaptive optics is another issue which may be consider in our case. 

a. Thermal Blooming 

Random temperature variations caused by turbulence take place when a 
high power optical beam travels through the atmosphere. As a result, the refractive index 
changes randomly and forms random lens along the beam path. This causes the optical 
beam to spread as it goes through the turbulent atmosphere. Thermo blooming is the 
laser-induced absorption-driven heating which spreads the beam core. Blooming is 
nonlinear and depends not only on medium composition and density but also on beam 
properties as well such as wavelength, pulse length and intensity. Adverse weather 
conditions, humidity, absence of cross winds and the smoke from combat environments 
can also degrade the beam quality. 

b. Absorption 

One of the major limitations of the transfer of laser energy from the source 
to the target is the absorption of the laser’s energy by the earth’s atmosphere. Light 
passing through an optical system can be attenuated by absorption as shown in Figure 4. 
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Figure 4. Light Absorption. 


The exponential law of absorption is the basic working relationship, but specific terms 
such as absorbance and transmittance are widely used as well. The differential absorption 
can be expressed as 

dd> = -ad>* (2.1) 

where a is the absorption coefficient and which upon integration from 0 to z gives the 
exponential law of absorption: 



d> 


0 


e 


-az 


( 2 . 2 ) 


c. Diffraction 

Diffraction limits affects propagation over a large distance and therefore 

has to be taken into account. For distances larger than the Rayleigh length, z > 7t / A,, 
where w is the radius of the source aperture, we are able to use the Fraunhofer 
approximation. At a wavelength of 1.0 |im and an aperture width of 2.5 cm (1 inch), this 
distance is z > 1250 meters. If we deal with a circular aperture, the intensity distribution 
can be written as [Ref. 13] 
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(2.3) 


J L kwr I z 

where A=kw^ , w is the radius of the source aperture, is the Bessel function of order 
1, and r the radius coordinate in the observation plane. The intensity distribution in 
Equation (2.3) is referred to as the Airy pattern after G. B. Airy who first derived it. Since 
sin0 = r! z, the irradiance can be written as a function of 0 , and Equation (2.3) gives the 

intensity 3 at distance z equal to 

3 = 43(0) 


/f (fcwsin0) 
(^sin0 f 


Eigure 5. 


The fundamental mode (TEMqo), from a circular aperture, is shown in 


3/3(0) 



- 7.5 -5 - 2.5 2.5 5 7.5 

/cwsinO 

Eigure 5. Intensity Pattern for Circular Aperture. 


To find the width of the central lobe, measured along the horizontal axis, 
the numerator of Equation (2.4) where it is zero (E' of minima) must be found. Erom 
tables, it can be seen that J^iu) = 0 when u = 3.83 so that the width of the central lobe at 
distance z is 

D = 1.22—. (2.5) 

w 
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Equations (2.3) and (2.5) strongly suggest short wavelength for a directed energy weapon 
as the intensity at the target is inverse proportionally of the square of the wavelength, 
while the width where the beam is focused is proportionally to the wavelength. 
Furthermore, circular apertures are preferable as 84% of light arrives within the central 
lobe (Airy disk), and 91% within the bounds of the second dark ring [Ref. 6]. For the 
FEE weapon with A- = 1 p-m wavelengths laser light from a w = 0.2 m aperture on the 
ship, the spot size at range z = 5 km could be as small as Wj = cm. To avoid thermal 
blooming, to damage a larger target area, the beam only will be focused to a 10 cm spot. 
Diffraction at A = 1 pm is not a problem. 

d. Scattering 

Rayleigh scattering refers to the scattering of light off molecules in the air, 
and from particles up to about a tenth of the wavelength of light. The strong wavelength 
dependence of Rayleigh scattering favors the short wavelength A since the scattered 
intensity 1 is proportional to A^ [Ref. 14]. For example the Rayleigh scattering at 400 nm 
is 40 times greater than at 1pm for equal incident intensity. 

For particles sizes larger than a wavelength, Mie scattering predominates. 
This scattering produces a pattern like an antenna lobe, with a sharper, more intense 
forward lobe for larger particles. Mie scattering is not strongly wavelength dependent. 
Both scattering types are presented schematically in Figure 6. 


COMBAT ENVIRONMENT 


Rayleigh Scattering MieScattering 



Mie Scattering, 
larger particles 

> 
-► 
> 
-► 



Incident light 


Figure 6. Scattering from Air Molecules. 
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The atmospheric transmittance over a wavelength range is a complex 
problem and extremely crucial to directed energy because it determines the choice of a 
suitable wavelength. Absorption-driven thermal blooming, scattering and turbulence can 
be major limitations to performance and can lead to the dispersion of the beam power. 

The atmospheric transmittance problem has been described by computer 
modeling codes such as FASCODE, LOWTRAN, MODRAN and MOLLY. All the 
above codes allow the user to insert weather conditions and to select from several 
different environments. MOLLY is a time-dependent computer simulation of adaptively 
compensated laser propagation through turbulence and thermal blooming. A graph of 
typical atmospheric transmittance versus wavelength using this simulation is shown in 
Ligure 7. Notice there are specific windows where the attenuation is much less, and 
therefore it is preferable to operate the laser in those regions. 


Maritime Atmosphere 
Marine Aerosoi, Horizontai Path 



Ligure 7. Atmospheric Transmittance (Lrom [Ref. 15]). 

As can be seen from the above figure, both COIL and MIRACL lasers 
perform poorly in a maritime environment and suffer from absorption. Wavelengths 
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A-=1.06|im, A-=1.25|im, A-=1.62|im are attractive for the MW-class because as the 
absorption coefficient is small, a < O.OOlkm'^ For the value of a =0.001 km"^ we have 
almost no reduction of the initial power, see Equation (2.2), and d>o . Note that the 
width of each window is relatively small (~0.2pm), and a small deviation causes 
excessive dissipation and spreading due to atmospheric conditions, which results in 
degraded range and beam quality. 

3. Adaptive Optics 

Current acquisition and pointing systems are capable of tracking maneuverable 
supersonic targets and holding the optical beam to a small area of-100cm. Adaptive 
optics can be used to transfer the lethal beam and focus it in. In addition, adaptive optics 
are used to overcome propagation affected by turbulence.. The advantage of using 
adaptive optics is shown in the simulations results in Figure 8. 



no adaptive optics adaptive optics 

Figure 8. Function of Adaptive Optics (from [Ref. 15]). 

On the left hand side, no adaptive optics are used and thermal distortions in the 
atmosphere spread the laser beam. On the right hand side the improvement is obvious 
after the application of adaptive optic mirror at source. Thermal distortions and moderate 
thermo blooming are largely corrected and allow the delivery of the lethal beam to the 
target. The system performance is degraded if the thermal distortions and the thermo 
blooming effect are more than moderate. In the FEE application discussed here, 
turbulence is much less than shown in Figure 7 so that adaptive optics may not be 
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required. If it is required the adaptive optics adds a cost of $3M to the entire FEL system 
[Ref. 8]. 

C. FEL SAFETY 

1. Radiation Concerns 

Radiation shielding is required for the electron beam dump. The FEL with the 
energy recovery system allows the electron beam to reenter the RE cavities 180 degrees 
out of phase with the accelerating fields. This energy is recovered as RE power and used 
to accelerate subsequent electron pulses to repeat the cycle. The residual energy after 
deceleration is about 7 MeV [Ref. 16] and is below the threshold of generating neutron 
radiation when it is dissipated in the beam dump. Only beryllium (Be) and Osmium (Os) 
have a neutron production threshold below 7 MeV and therefore, as long as beryllium 
and osmium are excluded from beam dump structural materials, neutron generation is not 
an issue [Ref. 17]. Without the production of neutrons, much less shielding is required. 
Other human safety issues include the need for eye protection when personnel may be 
topside or have visual access to the topside scattering off the sea surface during an 
engagement. In such cases, safety goggles may be required [Ref. 18]. 

2. Helium Concerns 

The cryogenic systems require the use of liquid helium. Helium is non¬ 
flammable, non-corrosive, non-toxic and not listed as a marine pollutant by Department 
of Transportation (DOT). When the FEL is running properly, the helium is in a closed 
loop refrigerator and totally contained. Thus, there are no safety concerns. However, 
because the helium is cryogenic, a vacuum is used to insulate the piping and storage 
containers (dewars). If a major vacuum leak should occur, the helium boils off quite 
rapidly leading to safety concerns: 

• It can over pressure the storage tank or piping since helium expands 750 
times when going from a room temperature liquid to a gas. If a suitable 
way is not provided for the gas to escape, the vessel could explode. 

• If the expanding gas should come into contact with someone, cryo-bums 
can occur because it is very cold when first released. 

• If the helium goes into a confined space, it can displace the oxygen and 
lead to suffocation of personnel in those compartments. If air with less 
than 12% oxygen is breathed, unconsciousness can occur in 9 seconds 
with brain damage occurring in 5 minutes. If the level drops to 8%, death 
is imminent in minutes. 
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The precautions are: 

• Provide burst disks on the system connected to a vent so the helium can 
safely escape. 

• In the case of an expanding plume, do not stick a body part into it. 
Protective clothing is required. 

• Provide a means for the helium to escape if it wants to rise. If there is a 
pipe in the ceiling, it will escape up the pipe and disperse. To prevent 
leakage into other compartments have a ledge at the top of the door below 
the ceiling level to prevent the flow. Also, oxygen monitors on the ceiling 
to send out an alarm if oxygen levels drop below 18% are required. 

The total volume of liquid helium in the linear accelerator will be less than 5000 
liters, making approximately 5 megaliters of gas available which could easily fill a room 
with 155m volume. Nevertheless, helium is much safer than liquid nitrogen. It goes to 
the ceiling and escapes. Nitrogen goes to the floor and cannot be seen once it warms [Ref. 
19]. Helium vents would be included in the weapon design. 

D. CONCLUSIONS 

1. Challenges 

Theoretical calculations have shown that FELs have the potential to scale to 
megawatts, but it is an engineering challenge to improve reliability and reduce the size 
and the cost of an EEL. Eor instance, a critical concern is improving electron-beam 
brightness as beam current is increased. Additionally, new functions must be introduced 
into EEE technology to advance the MW system. The radius of curvature of the 
resonators mirrors and a cooled sapphire window must be controlled to tight tolerances. 
Eurthermore, atmospheric propagation effects, and the thresholds at which these effects 
are significant, must be well understood. This understanding will, to a large extent, define 
the required laser weapon system. Performance, downsizing and cost are also crucial 
considerations and are illustrated in Eigure 9. 
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HEL Systems Considerations 



Figure 9. Acceptable System Diagram (After [Ref. 20]). 

After the workshop in Virginia in June 2001, the synopsis is that a MW-class 
infrared free-electron laser (IRFEL) appears to be a challenging, but feasible, proposition 
[Ref. 21]. 

2. Resources Required 

The DoD investment strategy is to fund those technologies that have the potential 
to penetrate and bypass identified technology barriers. In practice, in the face of funding 
pressures, inadequate funding is being provided to a wide variety of programs. These 
programs should be replaced with focused, sequential developments funded at the level 
of effort needed to create real progress. It is also believed that considerably more funding 
on the level of $100-150 million per year is needed [Ref. 22]. 

Currently, the money spent in this area in FYOl was $30 million, while the 
funding for the DoD high-energy laser programs was $474 million [Ref. 20]. Studies at 
TJNAF and simulations presented in Chapter IV have indicated that scaling the system to 
100 kW of power is feasible and straightforward. Funding is already being supplied by 
the Office of Naval Research (ONR) and additional funding is expected from the U.S. Air 
Force. 
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m. THE PHYSICS BEHIND EEL 


A. INTRODUCTION 

A free electron traveling undisturbed will not spontaneously emit energy. 
However, if the electron accelerates, it must absorb energy. If it decelerates, it must 
radiate energy. For example, if an electron increases its velocity in a particle accelerator, 
it does so by absorbing energy from an RF electromagnetic field supplied by the 
accelerator. Its kinetic energy is 

K^={y-m/ (3.1) 

where y is the Lorentz factor and me is the electron rest mass. The term is called rest 
energy (0.511 MeV). The Lorentz factory is given by 

Y= , , (3.2) 

Vl-P P 

where [3 =13 /c is the dimensionless electron velocity. 

In an FEL, the electron trajectory is directed with bending magnets into the 
undulator. Along the path of the undulator, the electron starts to wiggle in the transverse 
direction as soon as it passes through each alternating magnetic field. During the 
wiggling, the electron emits radiation into a narrow cone of opening angle "0 [Ref. 23], 
where 

i3=y‘ (3.3) 

This type of radiation should be distinguished from Bremsstrahlung or broad band 
synchrotron radiation which usually consists of several harmonics in addition to the 
fundamental [Ref. 24]. This thesis deals primarily with a pulsed FEL operated as an 
oscillator. A radio-frequency (RE) electron-gun is used to inject electron pulses into the 
accelerator. The radiation is stored in a cavity and amplified over many passes. Other 
types of EELs either use continuous wave (CW) power or amplify the radiation in a 
single pass and are called amplifiers [Ref. 25]. Table 1 lists not only the different types, 

but also the design parameters for existing and proposed EELs [Ref. 26]. 
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FA - FJectrostalic Accelerator 
O - Fl-F Oscillator 
S - SASF FFI. 


Table 2. The Short Wavelength PEL in 2000 (From [Ref. 26]). 
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The first column lists the operating optical wavelength X which is given by [Ref. 25] 


1 = 


2y" 


(3.4) 


where X,, is the undulator wavelength, and K is the undulator parameter listed in the 
following column. The relativistic Doppler shift and Lorentz contraction due to the 
electron motion is included in the calculation of the radiated optical wavelength in 
Equation (3.4). This is called the resonance condition and occurs when the electron slips 
behind the optical wave by one optical wavelength as the electron traverses an undulator 
wavelength. From the simple relationship of (3.4), it is seen that an EEL can easily be 
designed to operate at a broad span of wavelengths, including wavelength regions where 
there are no powerful light sources. For example, no conventional laser operates in 
ultraviolet or an extreme ultraviolet range of the electromagnetic spectrum (300 to 30 
nm) and no conventional laser operate in the far infrared (1mm to 20 pm) [Ref. 27]. 


The second column lists the normalized electron pulse length 



(3.5) 


where ^ is the electron pulse length and c the speed of light. From this column, it can 
clearly be seen that the majority of the FEEs create micropulses and only a few work CW 

The third column provides the electron beam energy E as it comes from the 
accelerator, ranging from a few MeVs up to 1 GeV in ELETTRA while the proposed 
electron beam energy is up to 30 GeV in the TESLA FEE facility. The next column shows 
the average peak current 1 ranging from a few milliamperes to 300A for the existing 
FEEs and up to 5000A for the proposed FEEs. 

The fifth and sixth columns are parameters related to the undulator design. The 
number of undulator periods is N, and the undulator wavelength is 



(3.6) 
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where L is the undulator length. Note that in a few cases, in the A^-column, there are 
multiple undulator sections, which are referred to in the FEL klystron undulator designs 
[Ref. 26]. The next column also depends on the design and is the undulator parameter 




eBXg 
2k nif^ 


(3.7) 


where e is the electron charge magnitude and B is the rms undulator field strength. Take 
special note that the values of ^ ~ 1. For ^ »1, there will be radiation in many higher 
harmonics [Ref. 24], and for K«l, FEL gain would be small. 

The last column lists the accelerator type. Note that the majority are RF linear 
accelerators and oscillators, as in the experiment and simulations described in the 
chapters that follow. 

B. FEL INTERACTION-PENDULUM EQUATION 

In this thesis, as well as in most cases. Coulomb forces between the electrons are 
small because the beam energy is large. Relativistic electrons streaming in the z direction 
and interact with a helical undulator magnetic field, 

B = B, (coskoZ,sinkoZ,0) (3.8) 


where B is the magnetic field of the undulator and ko = 27t / is the undulator 
wavenumber. The electron Lorentz force equation is [28] 

+ (3.9) 

dt m^c ^ 

Initially, we find the electron trajectory through the undulator by assuming that there is 
no radiation F = 0. The electron energy change is given by 

(3,10) 

dt m^c 


which means that dy ! dt is zero according Equation (3.10) for the F = 0 case , so that y 
would be constant. The term [3 x.S in Equation (3.9) is 

P xfi = x(-fio|3^ sinkoz) + K^o Pz cosk^z) + z(RoPx cosk^z) (3.11) 
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Then, Equation (3.9) can be separated into components and becomes 




dt jm^c 


sin Lz, 


(3.12) 


which by integration gives 


P. =-^%;-cos k,z 

ym^c 


(3.13) 


since dzldt = |3j,c, and where it is assumed that the ejection of electrons is perfect 
and the constant of integration is zero. Substitute =2ji I \ Equation (3.13) becomes 


Px = -T^^^^cos V 

zTty m^c 


(3.14) 


Define the undulator parameter as K = so that Equation (3.14) becomes 


|3, = -—cos^o^- 

Y 


(3.15) 


Similarly for the y-component, we have 


P„ =-—(sin^oZ). 


(3.16) 


The vector addition of the components in Equations (3.15) and (3.16) gives the electron 

—► 

velocity in the transverse direction = (P^,P ,0) . Thus, 


P^ =-(cos^oZ,sin^„z,0) 

Y 


(3.17) 


Now that we have the electron trajectory, we add radiation to find the microscopic 
motion. Taking into account the radiation created in the undulator, consider a circularly 
polarized plane wave for the optical field of the form 


= ^(cosvi/,-sin\|/,0). 


(3.18) 
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= £'(sin\|/,cosv|/,0). 


(3.19) 


Here V|/ = fe - co? + (p , £ is the electric and magnetic field amplitude in cgs units, k = 0) I c 
is the optical wavenumber, and cp is the optical phase. 

With radiation included Equation (3.9) becomes 

— =-^ECP^cosiiz-P sinv|/). (3.20) 

dt nif 

Substitution of Equations (3.15) and (3.16) into (3.20) yields 

pKF 

-cos(^+cp), (3.21) 

jm^c 


where y=dy/dt and ^+(p =\|/= (A:+^o)z-a)t+ (p . ^ is defined as an electron 
phase, which describes the electron position in the combined undulator and optical fields 
(^oand k terms). Note that in Equation (3.21), the rate of change in energy, y, is 
proportional to the electric field E and the undulator parameter K, but inversely 
proportional to the energy of the electrons y. As the energy of the incoming electrons is 
increased, the coupling with the optical field is reduced. 


Expanding Equation (3.2), the Eorentz factor reads 

y-^=l-P,^-Pl, 

and substituting Equation (3.17) in Equation (3.22) yields 

y-\\+K^) = \-^l. 

Taking the time derivative on both sides of Equation (3.23) gives 

- 2 y-iy(l+r-) = -2p,|3, = 

y \ + K 

Recall that the electron phase is 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


and by differentiating Equation (3.25) twice gives 
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^ =[(fc + fco)i]-CO = [(fc + ^o)Pz]-®’ 


(3.26) 


^ - (A:+A:o)(3^c ^ - 

{k+k^Jc 

and by inserting the result of P^. in Equation (3.24) we get 

Y_ Y'P.C 

Y i\+){k+k q)c 


(3.27) 


(3.28) 


For relativistic electron, k » A:,, and P^ ~ 1, so that Equation (3.28) can be written as 

t = (3.29) 

Y i\+K^)kc (l+iY")co 

The resonance condition of Equation (3.4), with % is replaced by 271 c/CO can be written 
as 


(0 = 


tOo2Y 
l + K^' 


(3.30) 


Substituting Equations (3.21) and (3.30) in Equation (3.29) and solving for ^ gives 

•• 2m pKF 

^ ^i^o£^cos(C+(p). (3.31) 

m^ci 

Dimensionless parameters are introduced to make the model simpler and more 
physically meaningful. The dimensionless time X =ctlL describes the electron and 
optical evolutions that take place as x ranges from 0^1. The time derivatives now are 

(..) = d(..)ldx instead of (..) = d(..)ldt. The dimensionless optical complex field is 
defined as 

a = \a\e“^ (3.32) 

II 2 2 

where |a| =4nNeKLEly me is the dimensionless optical field amplitude. Applying all 

the new dimensionless quantities. Equation (3.31) takes the form of a simple pendulum 
equation 
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(3.33) 


^ = |a|cos(^ +(p). 

Equation (3.33) governs the electron’s microscopic phase dynamics relative to the 
optical field inside the undulator. Therefore, the complicated motion of the electrons is 
described by this simple and well-known pendulum equation The positive or negative 
acceleration in Equation (3.33), is proportional to values of the cosine term and results in 
the bunching of the electrons within an optical wavelength. Consider a group of electrons 
uniformly distributed in phase. The electrons with phases in the range 
-7t/2<^ -rep <7t /2 will gain energy and speed up, while the electrons in the range 
7t/2<^-i-(p<37t / 2 will experience a loss in energy and slow down. As a result, the 
electrons have the tendency to bunch together near ^ =7t /2. This bunching coupled with 
the wave equation contributes to a net energy transfer to the optical field. 

The dimensionless field amplitude |a| of Equation (3.33) determines the electron 
bunching rate along the undulator. If |a|»7t, the optical field is strong and the bunching 
occurs quickly. If |a|«7t , the optical field is weak and the change in the electron phase 
velocity is small [Ref. 24]. 

The wave equation that drives the complex optical dimensionless field of 
Equation (3.32) can be put into the following simple form [Ref. 25] 

a = -j< e-‘^ > (3.34) 

where the dimensionless electron beam current j = SN(eKKL)^ p , p is the actual 

electron beam density, and (..) in the wave equation indicates an average over all the 
electrons. This average will be non-zero only when the electrons are bunched. The rate of 
change of the optical field is proportional to the dimensionless current j . 

The pendulum Equation (3.33) and the wave Equation (3.34) are coupled together 
through the dimensionless current y. If j < I, the gain is low, and if j »1 the gain is 
high. Both the pendulum and wave equations, even in this simple form, are valid in high 
and low gain regimes, as well as in weak and strong optical fields. 
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C. PHASE SPACE ANALYSIS 

The electron evolution and the optical wave are coupled accordingly in Equations 
(3.33) and (3.34). The phase-space evolution of sampled electrons is used to describe the 
motion of the electrons over an optical wavelength. Each electron is started with its initial 
conditions ^ (0) and =v (0) at the beginning of the undulator (x =0). The phase 
velocity is given in dimensionless notation as 

^L[(k+k,)^^-k]. (3.35) 

dl 

At resonance, v = 0, so solving (3.35) for (3^ and using the resonance condition 
X = \(l + K^) /2 y^ and the approximation that y »1: 



k 

k -hkQ 


l+K^ 


(3.36) 


Substituting (3 into Equation (3.35) yields 


V = L 


\+K^ 

ik+k,)i\—^)-k 

2y 


^-Lk 




(3.37) 


using kQ« k. Now differentiating Equation (3.31) with respect to y, substituting 
L = AA,) , and applying the resonance condition gain gives: 


Av =471A 



(3.38) 


Erom Equation (3.38), it is clear that the change in phase velocity is related to the 
change in the electron energy. Eigure 10 shows the phase space evolution for electrons at 
resonance and above resonance. In the first case, electrons are injected into the undulator 
with initial electron velocity v q = 0 . They appear as yellow spots (evolution denoted by 
yellow going to red). Note that the electrons between -7t/2<^ <7t /2 increase their 
phase velocity by absorbing energy from the optical field. An equal number of electrons 
between 7t/2<^ <3% / 2 decelerate and transfer energy to the optical beam. Thus, the 
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net energy transfer is zero. On the top right of the figure, where the output of the optical 
gain G(x) is located, it can be clearly seen that a resonant electron beam has no gain. 


**** FEL Phase Space Evolution **** 

j=l ao=7i: Vo=0 N=36 



(a) 


**** FEL Phase Space Evolution **** 
j=l ao=7r Vo=2.6 N=36 


Iz 


Gbin 


□ .3 



(b) 


Figure 10. 


The Phase-Space Evolution at Resonance (a), at Optimum Initial Velocity (b). 
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On the other hand, the bunching that occurs near the relative phase ^ +(p =7t /2, 
drives the optical phase cp in Equation (3.33). The optical phase (p(x) is plotted below 
the gain and it is relatively large, (~0.1). 

Electrons injected slightly above resonance at Vq = 2.6, are shown in the bottom 
part of Eigure 10. In the phase space evolution for x = 0 —> 1 along the undulator, 
bunching near ^ =7t occurs. This bunching indicates that a significantly larger number 
of electrons transfer energy to the optical field and the gain is increased to G = 12.5% . 

The curve is shown in red in the same figure is called the “separatrix”. It separates 
the closed and open orbits in phase space (phase-space paths) and connects unstable fixed 
points at ^ = -7t /2, 371/2 and V = 0 . The separatrix points are given by [Ref. 

25] 

vf = 2|a|[l+sin(C,+(p)]- (3-39) 

The peak-to-peak height of the separatrix is A\a\ . Electrons outside the 

separatrix follow open orbits; those inside the separatrix are trapped in closed orbits. 

D. GAIN 

Erom the wave Equation (3.34), it is clear that if 7=0, then a = 0 and there is no 
gain. If the dimensionless current is small j «Tl , then a small change in the optical field 

amplitude occurs, and the gain is low since the electron phase velocity does not change 
much during the evolution through the undulator. Today, the majority of the operating 
EEEs are low-gain, low-current EEEs, and based on the following assumptions, it is 
possible to derive the weak field gain. 

• Assume weak fields: a «%, cp = 0 

• Assume low gain: a-O^a-a^, (p~0 

With these assumptions, the pendulum equation can be solved analytically by expanding 
^ and V in a power series. 

The zero order expansion is 
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(3.40) 


C'°^=Co+VoX 

V® =Vo. 


To first order in , where a,, = |a(0)|, the pendulum equation gives 


^ ® =-^[cos (Co +VoX )-cosCo +VoX sinCo] 
V® =[sin(Co+VoX)-sinCo]- 


(3.41) 


In order to have net energy transfer from the electrons to the optical beam, we 
must satisfy (v)-Vo ^0, where (v) is the average phase velocity over all the electrons. 

The average of the in Equation (3.41) is zero because the average of the sin 


function is (sin) = 0 which means that just as many electrons gain energy as lose energy. 
Therefore, second orders terms are needed to obtain non- zero gain, so that we find 


V ® = -^ I (cos(2C 0 + 2voX) - cos(2Co)) + cos^' qX ) - 1 -v qX sin(C o)cos(C o +voX ) 


V 


1 


(3.42) 


The average of the second order term is 



=Vo+- 


0 L 


cos^ qX ) -1 -I- sin^ qX ) 


(3.43) 


The above dynamics of the electron position is related to the energy transfer, with the 
relation Av =4'KN{Ay /y). Both sides are averaged to find the average change in the 
electron beam energy 


Ay mc^ =ymc‘ 


((v)-v.) 

471 A/' 


(3.44) 


Energy conservation is used to calculate the gain G(x). Gain is defined as the ratio of the 
energy transfer from electrons. Equation (3.44), divided by radiation beam energy in cgs 
units, 2E^dV I Sk, in a volume dV, where is the electric optical field. The filling 
factor F is defined as the ratio of the electron beam area divided by the optical mode area 
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(3.45) 


where is the electron beam radius and Wg is the optical beam radius at the waist. The 
number of electrons in a volume dV is given by pFdV , where p is the electron density. 
Thus the gain at time x is given by 

ioFdV )ymc^{{v)-vAl AnN 

G{x)=-— - - —-. (3.46) 

2E^dV I %K 

In the exchange of the energy, the second order contribution of Equation (3.43) is used. 
Also, the definition of dimensionless current j = SN(eilKL)^ p Fly^mc^ is used. Using 
these relations in to Equation (3.46) yields 


G(x) =^{2-2cos(VoX)-VoX sin^gX)}. (3.47) 

Vo 

Erom the above equation, the strong relationship between the gain of the optical field and 
the initial phase velocity Vg can be seen. 


Eigure 11 shows the output of simulations in the weak-field, weak-current regime. 
The final gain G(x = 1), as well as the final optical phase shift (p(x =1) versus the initial 

phase velocity, Vg, is plotted. 
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Figure 11. Weak-Field, Low-Current Gain G(V(,)and Optieal Phase Shift cp^o) Speetra. 

The simulation result is in total agreement with the analytie results of Equation 
(3.47). The initial phase veloeity, whieh is a funetion of the eleetron beam energy (recall 
that Ay /y = AVq ! ‘XtzN), determines the final gain. At the resonance condition, (Vg = O), 

the gain is reduced to zero. On the other hand, the optical field drives the optical phase 
shift, with a peak value of Acp = O.Sy . There is a peak gain of G = 0.13 j when the initial 

velocity is Vg = 2.6 and the optical phase shift is reduced to Acp =0.02. As the gain 
spectrum is anti-symmetric around Vg=0, there are losses for negative initial phase 
velocities (Vg <0) and accordingly, the maximum absorption occurs at Vg =-2.6 . 

E. SHORT PULSES 

From the last column of Table 1 it is clear that the majority of FELs use radio¬ 
frequency accelerators (RE). Electrons are ejected in short pulses from an electron-gun 
with a frequency / and are accelerated. The separation distance between the pulses is c/f. 
Thus, for /= 750 MHz, the pulses are separated by 40 cm. The pulse length that is 
currently used for FEE oscillators is often picoseconds and is comparable to the slippage 

distance NX, where N is the number of undulator periods. 
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In the picosecond scale, the short-pulses effect the FEL interaction. As electron 
“buckets” enter an undulator length L , short optical pulses are produced due to 
spontaneous emission. These optical pulses bounce between the mirrors of the resonator 
separated by distance S. The optical pulse after a time 25 /c arrives back at the 
beginning of the undulator. If the incoming electron buckets enter the undulator at the 
same time and coincide with the optical pulse, exact synchronism occurs. The 
displacement between optical and electron pulses is called desynchronism, d, and in exact 
synchronism d = 0. As the optical pulses enter the undulator simultaneously with 
electron pulses,gain development is delayed. This effect is called lethargy [Ref. 25]. The 
electron pulses travel slightly slower than the optical pulses. The electron pulses fall 
behind and preferentially amplify the trailing part of the optical pulse; so the centroid of 
the optical pulse travels slower than the speed of light. After each pass, the optical pulse 
falls farther behind the electron pulse and consequently, after some passes, the power 
starts to decay. 

Figure 12 shows the short pulse evolution at exact synchronism, d = 0 . The lower 
left window shows the electron pulse position for x = 0 and X = 1 . The electron pulse 
shape is parabolic and has the form j(z)- j{l- 2z^ for j(z)>0 and is zero 

otherwise. The normalized pulse length is o^and the dimensionless current j = 5. The 
middle-left window corresponds to the field evolution, |a(z,n)|, over the passes n=500 
passes. Note that the optical field shifts to the left as the centroid of the light pulse travels 
slower than the speed of light c. The final optical pulse shape |a(z)| is shown in the upper 
left window. In the specific simulation output, the optical field and the power decay 
accordingly after n ~ 150 passes. 
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Figure 12. Short Pulse Evolution at Exaet Synehronism. 

The middle lower window represents the weak field spectrum G^„)and is plotted for 
reference. The window above that shows the evolution of the optical power spectrum 
P(v,n)over n = 500 passes. The final power spectrum P(y) is in the top-middle 
window. The pointed tick-mark at the top indicates the initial optical wavelength derived 
from the resonance condition. The thicker mark indicates the center of the final spectrum. 

The lower-right window shows the dimensionless optical power P{n) at the end 
of each pass. The middle right window shows the evolution of each of the electron energy 
spectrum f(y,n) over n = 500 passes. The upper-right window shows the final electron 

spectrum after n = 500 passes. Again, the pointed tick-mark shows the initial phase 
velocity of the electron beam at resonance. The thicker tick-mark indicates the position of 

the final average electron phase velocity (Av ). 
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In exact synchronism, the steady state power is reduced to zero. The sbwer speed 
of the light pulse centroid is overcome by reducing the path inside the resonator by A5 . 
To do so, piezoelectric crystals move the mirrors slightly and d =-2AS / NX. Figure 13 
shows, an FEL has the same design characteristics as the FEL in Figure 12, but with 
d = 0.06 instead of d = Q . 


7=5 0^=1 j=o.06 

2=10 ^ = 36 5^=0.0001 



-1.5 z 1.5-47 V 47 0 n 500 


Figure 13. Short Pulse Evolution at Relative Large Desynchronism. 

The introduction of a small amount of desynchronism d = 0.06 allows the FEL to 
produce a steady state power after n = 50 passes. 

As will be seen in the chapters to follow, desynchronism is important to the 
characteristics of pulse evolution. By increasing desynchronism d, less power is obtained 
but the FEL becomes more stable. The final optical pulse is longer than the electron pulse 
and as a result, the power spectrum is narrow. Increasing the desynchronism too much 
will cause the electro-optical pulses to not overlap sufficiently and consequently, the FEL 
will not operate. 
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IV. SIMULATIONS OF THE lOOKW TJNAF EEL 


A. INTRODUCTION 

The Thomas Jefferson National Accelerator Facility (TJNAF) FEL has 
demonstrated 2 kW output power of tunable infrared radiation [Ref. 29]. This is the 
highest average power level from an FEL in the world and was achieved with a design 
leading towards the further development of get higher power. A schematic representation 
of the system is shown in Eigure 14. 



In order to operate at higher output power levels, the TJNAE EEL utilizes electron 
beam recirculation. Energy recirculation or “energy recovery” is the process by which the 
same electron beam is returned to the accelerator. After the interaction inside the wiggler, 
the recirculated electrons arrive back at the linear accelerator “Linac” 180 degrees out of 
phase so that their energy is converted back into RE power used to accelerate subsequent 
electrons. 

The energy spread that a specific system can recirculate is limited. The limitations 
come from the induced RE phase shift upon arrival at the wiggler, which is equal to 
optical desynchronism. The change in desynchronism alters the FEL gain and the laser 
output power. Consequently, these optical power changes result in a change in the energy 
of the recirculating electron beam, which can potentially lead to additional electron beam 
loss on apertures and a phase shift of the decelerated beam. The power variations, the 
phase shifts and the beam loss can change the beam-induced voltage in the accelarator 
cavities. The process can be refered to as a “beam loading instability” [Ref. 31]. If the rf 
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control system does not possess sufficient gain and bandwidth, the whole recirculating 
PEL becomes unstable. 

Another constraint on recirculation occurs in bending magnets of static ble 
magnetic fields are used for the 180 degrees turns. The momentum of the electrons 
spread resulting from the EEL interaction leads to a range of bending radii that can drive 
some of the electron beam into the wall of the vacuum pipes and cause the EEE to cease 
lasing. Overall , there are many benefits to recirculation, such as the enhancement to the 
wall-plug efficiency and reduction of harmful radiation in the beam dump, so 
recirculation is desirable. 

In the course of developing megawatt power levels, which is required for a 
weapon, it is advisable to proceed step by step. Studies at TJNAE have indicated that the 
system can be upgraded to operate at an increased power of 100 kW. This upgrade has 
been proposed with a projected demonstration scheduled for 2005 [Ref. 11]. 

To support this goal, simulations were conducted to determine gain, steady-state 
extraction efficiency and the induced energy spread. The input data for simulations were 
design parameters chosen by the TJNAE scientists [Ref. 32]. To achieve 100 kW average 
power, the kinetic energy of the electron beam increased to = 210 MeV and the pulse 
repetition rate to = 750 MHz. The average electron beam power is given by 

= (4.1) 

where 7 = 270A is the peak current in the electron micropulse, 7^ =0.1 mm is the 

electron micropulse length and c is the speed of light. The average electron beam power 
then calculated as P^ = 14MW. Therefore an output power of 100 kW requires an 
extraction efficiency of approximately T] = 0.7 % . 

The undulator period is Xq = 8 cm over N =36 periods with an rms undulator 
parameter of K = 1.1. Table 1 summarizes the experimental parameters. 
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Electron energy 

£’,= 210MeV 

Peak electron pulse current 

/ = 270 Amp 

Electron pulse length 

fg = 0.1 mm 

Number of undulator periods 

A = 36 

Undulator wavelength 

= 8 cm 

Undulator length 

L = 288 cm 

Resonator length 

5 = 36 m 

Optical wavelength 

A- = 1 mm 

Step-Taper 

A = 0,±7t,±27t 


Table 3. TJNAF lOOkW FEL Parameters. 


In the simulations presented here, different undulator designs were investigated to 
aeeess their possible advantageous effects on the FEL interactions. The conventional 
undulator possesses a periodic magnetic field and wavelength, but a linearly tapered 
undulator gradually changes the undulator parameter K as shown in Eigure 15 by 
modifying the gap between the undulator magnets along the undulator length. The field 
can either increase or decrease along the undulator length. 


LINEAR TAPER, 6 
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Figure 15. 


Schematic for a Linear Taper Wiggler. 
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A step-tapered undulator abruptly changes the value of the field and K halfway 
through the undulator as Figure 16 shows. The stepped field can decrease or increase as 
shown. 
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Figure 16. 


Schematic for a Step-Taper Wiggler. 


The resulting PEL interaction in a tapered undulator [Ref. 25], [Ref. 33]-[Ref. 37] 
is described by the modified pendulum equation 


v=^’=5 +0 


X — 

V 2, 


A-l-|a|cos(^ -l-tp) 


(4.2) 


[ 0 for z < 0 

where 0(z) = ( is the step function, 

1 for z > 0 


={k+k^)z-(iit is the electron phase, k=2%lX is the optical wavenumber, 
k^ =2ji is the undulator wavenumber, ()) = kc is the optical frequency, V is the 
electron phase velocity, cp is the optical phase. 


I I ^'9 
a =\a\e 


(4.3) 


a 


is the complex dimensionless optical field, |a| = 2%NeKL\E\/j^mc ^, and. 


d^-4nNK^ 




j 


(i+i^ ) 


(4.4) 
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where 5 is the phase acceleration caused by the linear tapering of the undulator, 5K/K 
is the fraction change in the undulator parameter along the whole undulator length, 


A^-4nNK^ 


AK ^ 


(1+^) 


is the step-taper along the undulator at 


X = 


z 

NK 


2 ’ 


(4.5) 


(4.6) 


and AK IK is the step change in the undulator parameter at x =1/2. 

B. WEAK FIELD GAIN 

The 100 kW TJNAF FEL is described by the dimensionless current j = 5, the 
electron pulse length 


= 



(4.7) 


and the resonator cavity quality factor of Q= 4.2 . To study the FEL gain, a weak optical 
field is used so that the dimensionless field is |a| < % The desynchronism d is given by 


d = 


A5 

NX 


(4.8) 


and measures the shortening of the resonator cavity length by AS compared to the 
slippage distance NX. The desynchronism between the optical pulse and the electron 
pulse was varied from d = 0 to d =0.3. The EEL gain results from many simulations at 
different values of desynchronism d with a step-tapered undulator are plotted in Eigure 
17. 
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Figure 17. Weak Field Gain G versus Desynchronism d. 


For the conventional non-tapered case, (A = 0), it can be seen that the highest gain 
is 60%. At small and large values of desynchronism, the gain decreases for all the 
undulators away from peak values around d =0.1. The gain for small step-tapers of 
A = ±71 is only slightly reduced from the conventional case. Larger step-tapers A = ±271 
cause a significant reduction in gain down to around 40%, peaking at the lower values of 

d -0.05 . Also, the FEL still works above the threshold losses {Q=4.2) for larger 
values of d < 0.2. 

Figure 18 shows the gain spectrum for an FEL with an inverse taper rate of 
A = -271 over N = 36 periods with a current density of j = 6 . These parameters are 

descriptive of the proposed TJNAF 100 kW FEL. 
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Figure 18. Gain vs. for Inverse Step Taper A = -27t. 


For a negative step taper, the gain decreases in strong fields but still remains 
significant. Also, increasing the optical field ao results in a small increase in the phase 
velocity for peak gain. 


In the electron pendulum phase-space, determined by Equation (4.2), the negative 
step taper A = -2jr FEE shows 1% efficiency with 3.8% induced energy spread. Eigure 
19 shows a simulation in a strong optical field where electron bunching is evident in 
phase space evolution. 


**** FEL Phase Space Evolution **** 
j=5 ao=15 vo=7.4 N=36 

00=1 A=-2ji 

12| 0*iii 0.25 



0.02 



Eigure 19. Phase Space for Step-Taper A = -27C Undulator. 


45 







As can be seen in the figure the electrons start at initial phase velocity Vg =7.4 
with a random spread in phase velocities of width = 1. The light colored dots show the 
phase space position at the beginning of the undulator, while the darker colored dots 
show the position at the end. At phase position ^ -iz, the electrons are bunched, amplify 
the optical field and lead in a gain G -23%. The optical phase at this position is small 
(|) -0.02 and further reduced to (|) -0 at the end of the undulator. 

C. STEADY STATE EXTRACTION EFFICIENCY 

For a conventional undulator with no taper, the desynchronism was varied from 
d =0.005 up to 0.3 to study the pulse slippage effects. The best results for 
desynchronism were d =0.04, with an efficiency of T] =0.08 above the requirement for 
lOOkW. At that value of desynchronism, the induced energy spread was 
Ay /y = Av /4 k N = 3%, which is well below the 15% limit for recirculation. Figure 20 

shows that slight increases in desynchronism to d > 0.06 makes the efficiency marginal 
(r| < 0.7 %) for reaching 100 kW. 



Figure 20. Efficiency p Versus Desynchronism d. 
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Additionally for a step-taper undulator, the results of longitudinal multimode 
simulations are shown in the same figure. Again, the efficiency is plotted versus 
desynchronism d for different values of A . For the same value of d =0.04, a negative 
step-taper A = -71 once again demonstrated the best results, and the efficiency increased 
to-0.9%. 

Figure 21 presents a high power example for the A = -71 case with a 
desynchronism value of d = 0.04 . 

H<H<H<H<H<H<H< FEL PulSG EVOlUtiOH H<H<H<H<H<H<H< 

7 = 5 a, =3 c/=0.04 

g = 4.2 A = -71 A = 36 



-6 z 6 -50 V 50 -50 V 50 



-6 z 6 -50 V 50 0 n 500 


Figure 21. Power Simulation for Inverse Step-taper Undulator. 

The efficiency for this case is T] =0.86% , which corresponds to an output power of 116 
kW. The induced electron energy spread is only Ay /y = 2.4% , which is within the 15% 

specification for recirculation. The power reaches a steady state in less than n = 50 
passes and remains steady throughout the whole experiment. The good behavior of the 
power evolution P(n) is indicated at the bottom right window (steady state power after 
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500 passes). Larger values of tapering A<-27U, and desynchronism d >0.06 did not 
work. 

With both the linear and step taper, a slight increase in efficiency over the 
conventional undulator was observed for the small negative taper. 

The induced energy spread did not change much with any undulator design as is 
shown in Figure 22. 



desynchronism d 


Figure 22. Energy Spread Ay/y versus Desynchronism d for Step-Taper with Q=4.2. 
Lower Power Reduces Energy Spread for all Undulator Designs. 

The maximum values of the induced electron energy spread were found in the region 
with maximum efficiency. Taper rates A = 0,-7t result in energy spreads of 

Ay/y =3.5% and 3.7% respectively. These values were the largest found for all 
undulator designs and were far less than the maximum allowable limit for safe beam 
recirculation around 15%. Eor this reason, the quality factor Q was increased from 4.2 to 
10. The new resonator cavity will only have 10% losses per pass. Also, larger efficiency 
and induced electron energy spreads were examined. 
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D. HIGHER Q EXPERIMENTS 

The purpose of these simulations was to explore regions beyond the design 
parameters and find the optimum undulator design. The quality factor Q was increased 
from 4.2 to 10 while the rest of the input data remained similar. Such experiments have 
been conducted at TJNAF and are described in the following chapter. 

I. Step-Taper Wiggler 

One of the proposed designs was to make use of the step-tapered wiggler. For this 
design the extraction efficiency and the induced electron energy spread were explored. 
a. Extraction Efficiency 

In the simulations, the step-taper rates |A| were applied from 0 —> 471. In 
Figure 23 the curves show the result of simulations with a quality factor of 2 = 10. 



Figure 23. Efficiency p versus Desynchronism d for Step-Taper, and Higher Q=10. 
Power Above 100 kW for a Larger Range of Tapers A = 0,±7i,±27i. 


In this case, the extraction efficiency for some values of A and d exceeds 
1%. In this case of A = -71 , the efficiency increased from 0.9% to 1.4%. For all values of 
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A = ±471 and for a desynchronism of d >1.2, the FEL failed to exceed the efficiency 
limit. Compared to the case of Q = 4.2, there is a larger range of desynchronism in which 

the FEL operates above the extraction efficiency requirement of T] >0.7% . Moreover, 
other taper rates, such as A = +7t,±27t , were able to produce output power of exceeding 
100 kW. This result is significant, because it was not confined to a small region of 
parameters. See Figure 7 for 2 = 4.2. More acceptable taper rates of A and a larger d 
lead to more stable, and therefore, more desirable systems. The simulation output for the 
optimum inverse taper of A = -7U at d = 0.03 is presented in Figure 24. 


******* 

FEL Pulse Evolution ******* 

7=5 

Q 

II 

d=0.03 

Q 

II 

o 

> 

II 

A = 36 



—6 z 6 -50 V 50 0 n 500 


Figure 24. Power Simulation for A = -n,Q =10 . 

In this figure, the steady state power is evident and the induced energy spread is limited. 
The efficiency t] = 13.6% results in an output power of 190 kW which is almost double 
that of the objective. Also, for this power, the induced energy spread Ay /y = 4.3% is 
again within a safe limit. 
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b. Induced Energy Spread 

The new design in the eavity 2=10 inereased the induced electron 

energy spread significantly. Nevertheless, as Figure 25 shows, it is still below the 
maximum acceptable limit of 15% and our design can recirculate the electron beam. 



Figure 25. Energy Spread Ay/y for Step-Taper, 2=10. 

2. Linear Taper Wiggler 

A linearly tapered wiggler gradually changes the undulator parameter K by 
modifying the gap between the undulator magnets. A schematic representation was given 
in Figure 15. 

a. Extraction Efficiency 

In simulations, the linear taper rate 6 applied was varied over the values 
5 = 0,±7t,±27t,±47t, and corresponds to AK/ilf =±0.9%,±1.9%,±3.8% . Recall from 

equation (4.4) that 5 = -47t NK^(AK / .S:)(l+ ). 

The operating rate of desynchronism d between the optical and electron 
pulses was varied from 0 —>0.3. Figure 26 shows the extraction efficiency as a function 
of desynchronism and the taper phase acceleration 5 . 
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Figure 26. Efficiency r] versus Desynchronism d for Linear Taper with Q=10. Power 

Above 100 kW for a Larger Range of Tapers. 

Once again, it is found that the negative taper results in maximum 
efficiency. For the linear taper rate 5 = -2% , there was a maximum efficiency of 
T]-1.4% at <i=0.01. For the untapered undulator 6=0, there was a maximum 
efficiency of T] = 1.2% at the same value of d. Other undulators with negative taper rates 
of 5 =-7t ,-4k reach r] =1.3%, which is better than the zero tapered undulator. Note 
that 1.0% efficiency is easily achieved from almost all taper rate values of 6 with the 
exception of positive tapers 5 = 27t,4;t . For these taper values, r\ =0.9% and r\ =0.8% 

occurred respectively. These values of extraction efficiency are the smallest, but still 
result in an output average power of more than 100 kW. 

b. Induced Energy Spread 

A smaller exhaust energy spread allows the intense electron beam to 
recirculate without damaging components. Figure 27 shows the final induced energy 
spread produced by seven different designs. 
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desynchronism d 


Figure 27. Energy Spread Ay /y for Linear Taper, 2=10. 

The inverse taper undulator for 6 = -271 has a peak energy spread of 
Ay /y = 6.3% occurring at d =0.01. All the values of all taper rates, including the latter 
value, produce an energy spread below the maximum tolerable limit. An increase in 
desynchronism in order to reduce the electron beam energy spread is not necessary. Note 
that a positive linear taper produces the smallest energy spread compared to all other 
designs. In small values of desynchronism, all the taper rate values are almost 
comparable. However, after d = 0.04, the values that occur from the positive taper are 
significantly smaller. 

E. CONCLUSIONS 

The best results of 100 kW simulations for various designs and a quality factor of 
2 = 4.2 are summarized below in Table 4. 
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e = 4.2, STEP TAPER DESIGN 

Step-Taper 

Desynchronism 

Efficiency 

Induced energy 

Output Power 

A 

d 

ri(%) 

spread Ay/y (%) 

P (kW) 

0 

0.04 

0.82 

3.5 

115 

71 

0.04 

0.55 

3.2 

77 

271 

0.05 

0.25 

2.4 

35 

-71 

0.04 

0.86 

~3.5 

120 

-27C 

0.05 

0.64 

3.2 

90 


Table 4. Results of TJNAF 100 kW Step-Taper Undulator, for Q = 4.2. 

Step-taper wigglers were explored and compared with conventional wigglers. For 
the given TJNAF lOOkW FEL parameters, the inverse taper undulator A = -Tt achieved 
the highest final power at a desynchronism of d = 0.04 . The output power produced was 
F~120 kW which is above the 100 kW objective. The conventional untapered 
undulator, (6 =0), produced a final power of P = 115 kW, which again was above the 
objective at the same desynchronism value. In both cases, the induced electron energy 
spread Ay /y ~ 3.5% was smaller than the maximum allowable limit of 15% and did not 
limit the operational range of the FEL. 

The small values of Ay /y allowed for the possibility to explore a design with 
10% losses (2 =10) in an optical cavity. The high 2 =10 factor was chosen in such a 
way as to increase the efficiency of the system, while maintaining the induced energy 
spread well below the 15% engineering limit. The results related with Q, this design and 
the step-taper design are summarized in Table 5. 
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Q =10, STEP TAPER DESIGN 


Step-Taper 

Desynchronism 

Efficiency 

Induced energy 

Output Power 

A 

d 

ri(%) 

spread Ajlji^o) 

P (kW) 

0 

0.01 

1.19 

5.3 

167 

71 

0.01 

0.94 

5.1 

132 

271 

0.01 

0.78 

5.6 

109 

47t 

0.02 

0.12 

1.5 

17 

-71 

0.01 

1.36 

5.1 

190 

-27C 

0.01 

1.25 

5.3 

175 

-47C 

0.03 

0.45 

3.5 

63 


Tables. Results for Step-Taper, 2=10 Design. 

The best result for all the above eases was a taper rate of A = -7t and occurred at 
desynchronism d = 0.01 with an output power of P = 190 kW. Moreover, all the designs, 

except A = ±471 , operated above the threshold and are capable of producing an output 
power of at least 100 kW. The linear taper results are summarized in Table 6. 


2=10, EINEAR TAPER DESIGN 


Einear Taper 

Desynchronism 

Efficiency 

Induced energy 

Output Power 

5 

d 

ri(%) 

spread Ajlj(%) 

P (kW) 

0 

0.01 

1.19 

5.3 

167 

K 

0.01 

1.1 

5.5 

154 

2n 

0.01 

0.92 

5.1 

129 

4k 

0.01 

0.68 

5.5 

95 

-71 

0.01 

1.31 

5 

183 

-27C 

0.01 

1.43 

5.2 

200 

-471 

0.01 

1.32 

5.1 

185 


Tabled. Einear Taper, 2=10 Design. 
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The linear taper seems to be a very promising design. All taper rates, except for 6 = 471 , 
achieved 100 kW. For 6 =-271 at a desynchronism of d =0.01, the output power 
P = 200 kW was double the objective. From the table, the superiority of the inverse taper 
is clear. This kind of taper results in a large efficiency while the induced electron spread 
is comparable to other designs, and sufficiently lower than the 15% limit for safe 
recirculation. 

The common behavior for all forms of taper rates is a reduction in extraction 
efficiency and the energy spread for a slight increase in desynchronism. 

This work was presented at the 23rd International FEL Conference in Darmstadt, 
Germany and has been accepted for publication in Nuclear Instruments and Methods in 
Physics Research. 
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V. THE TJNAF TAPERED WIGGLER EXPERIMENT 


A. INTRODUCTION 

Using the Thomas Jefferson National Accelerator Facility (TJNAF) free-electron 
laser, experiments have been conducted which study its operation as a function of various 
taper rates and desynchronism values [Ref. 38]. The kinetic energy of the electron beam 
was = 34.5 MeV, the optical wavelength was A- = 6 |im, the electron pulse duration 

was 0.5 ps (pulse length =150|im), and the peak current was I -50A. The purpose 
of the simulations presented here is to compare experiment with theory, and, by using a 
wider range of parameters than allowed by experiment, to extract more physical meaning 
from the results. 

B. PARAMETERS AND SIMUUATIONS METHODS 

The simulations used the same parameters as TJNAF, described by dimensionless 
current density y = 10, normalized electron pulse length a^ = NA = l and undulator 

parameter K = 0.98(rms). Table 7 summarizes the experimental parameters and Table 8 
the corresponding dimensionless parameters that are used in the simulations. The 
desynchronism between the optical pulse and the electron pulse was varied from d = 0 
up to the highest value of d which produced power corresponding to an efficiency p > 
0.01%. As an example, the maximum desynchronism for the conventional non-tapered 
case (5=0) is d = 0.4. Not only were taper values for which we had experimental 
results explored, i.e. 5 = 0, lOrt, -Vtt , but values well beyond these were also explored. 

By using a wider range of values of linear taper, it was possible to plot the operating 
width of the desynchronism curve versus linear taper. The efficiency p and the induced 
energy spread Ay/y are determined as functions of desynchronism d and taper 
acceleration 5 . 
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TJNAF EXPERIMENTAL PARAMETERS 


Electron Energy 

Ymc^=34.5 Mev 

Undulator Length 

L= 1.1 m 

Undulator Periods 

N= 41 

Undulator Wavelength 

\ =2.7 cm 

Peak Undulator Magnetic Eield 

fi=5.5kG 

Peak Current 

7 = 50 A 

Taper Phase Acceleration 

(5/;t)=-7,0, 10 

Cavity Losses 

2=50 


Table 7. Parameters used in the Experiment 


DIMENSIONLESS SIMULATION PARAMETERS 

Current density 

7=10 

Electron Pulse width 


Undulator Parameter 

K= 0.98 

Desynchronism 

d= 0^0.4 

Taper Rate 

(5/7t) = 0,±3,±6, ±7,±8, ±9,±10,±12,±14,±16,±18,±20 +22,±24 


Table 8. Parameters used in the Simulations 
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C. UNTAPERED UNDULATOR 

For an untapered undulator, the desynchronism was varied from d = 0.002 up to 
0.4 in order to study pulse slippage effects. It was found that for desynchronism 
d = 0.02, the efficiency was T] = 0.04. However, at this value the power P(n) is not 
steady, as shown in Figure 28 and the power oscillations are caused by the trapped- 
particle instability resulting from strong optical fields [Ref. 25]. 

**** FEL Pulse Evolution **** 

j = 10 Oz=l d = 0.02 

Q = 50 N = 41 5^=0.0001 



Eigure 28. Trapped Particle Instability. 


In the upper left window the final optical pulse shape after n = 500 passes is 
shown. There, short optical pulses, modulated with sharp spikes, create two sidebands in 
the optical power spectrum P(v,n) drown in the upper-middle plot. In the upper right 
window, the electron spectrum f(v,n) shows the large induced spread in phase velocities 
Av =44, which is not desirable for recirculation of the electron beam. The 
corresponding induced energy spread 


Ay _ Av 
y AnN 


44 
471A 


= 9% 


(5.1) 
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Steady-state power without the trapped-particle instability or limit-cycle behavior 
[Ref. 39] is shown in Figure 29 for larger desynchronism ofd >0.15. Note the narrow 
electron distribution and the absence of sidebands in the power. Also note that the power 
evolution P(n) (bottom right) has no oscillations. 


**** FEL Pulse Evolution **** 
j = 10 az=l d = 0.15 
Q = 50 N = 41 8^ = 0.0001 



-6 z 6 -47 V 47 0 n 500 


Figure 29. Steady State Power for Conventional Undulator. 


D. GAIN SPECTRUM 

The dimensionless electron phase velocity V . is the meeting point between 
simulations and experiments and in dimensionless notation is given by 


V = 27tN 


(l + K") 
2y'X 




(5.2) 


where N= 41 is the number of undulator periods and y is the relativistic Lorentz factor. 
Single mode simulations are used to find the optimum initial phase velocity for 
maximum gain. 
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Figure 30 presents a plot of as a funetion of the tapering rate 5 . The relation 
for peak gain v„ = 2.4-5/2 is presented with the red circles and is valid in weak fields 
and moderate taper. Taper rates 5 >1071 and strong fields \a^\ >% change this relation. 



Figure 30. Plot of Phase Velocity versus Tapering Rate. 

Figure 31 presents the peak gain as a function of taper by using the optimum 
initial phase velocity of Figure 5. Is required for all tapering rates -247t <5 < 247t, 

we have gain G that exceeds the resonator losses. The optical cavity losses, used in the 
experiment, were for Q = 50 so that, gain G > 2% is required for operation. The 

minimum gain, after introducing taper was G = l%, and was met for 5 = 2471 . 
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The G(5) curve is almost symmetric in 5 with peak value at 6 = 0 where the gain 
is G = 160%. As the magnitude of 5 increases, the peak gain decreases and for 
5 = ±2471 reduces to G ~ 8% . There is only a 1% difference in gain at the ends of the 
curve where negative taper results in slightly more gain. 

Figure 32 (a) shows single mode gain plotted as a function of the initial optical 
field ao and the initial electron phase velocity v„. For 6=0, weak fields (ao< 7l) and 

current 7=10, the maximum gain is G~ 160% and occurs at v = 2.4. Figure 32 (b) 
shows that the application of modest positive taper rate (6 >0) makes the peak gain 
decrease and shift to the left in v„. Figure 32 (c) shows similar results for the inverse 
taper 6 <0but shifts to a reduced peak gain to the right in v^. For both positive and 
inverse tapers, it can be seen that gain is not symmetric around the resonance (Vq = O) and 

for large values of optical field strength ao, there is a distortion in the shape of the curve. 
Mention that positive taper has one gain peak in strong fields while negative taper has 
two, nearly equal, gain peaks in strong fields. 
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Figure 32. 



ln(l+G) 



j=10 N=41 


0.44 


▲ 

0.43 

j=10 N=41 

6=-8x 0=0 * 


(b) 


(c) 

G(Vo,ao) for Zero taper (a), 5 =871 (b), 5 = -Sti (c). 
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E. STEADY STATE POWER 

In our multi-mode simulations, the following values of linear taper were used: 

5 = ±3Jt +6k +7% ,±87t ,±97t ,±10rt ,±127t ,±147t ,±167t ,±187t ,±207t ,±227t ,±247t. 

Figure 33 summarizes the results of multimode simulations and plots the 
efficiency p for positive taper. The untapered undulator reaches its maximum 
efficiencyp =4% at d =0.002. As the taper rate increases the efficiency decreases. For 
5 =97U, the efficiency is 50% smaller and peaks at 0.006. For 5 =14;t, the 
efficiency diminishes significantly to p =0.8% at d =0.006 . After 5 =147t, the peaks in 
the efficiency curves stay at lower values of desynchronism. Sidebands are formed due to 
the trapped-particle instability and the efficiency is enhanced up to 2% for 6 = 1671 at 
d =0.002. Finally for larger taper rates 5 = 247t,the efficiency drops to only 1% at 

J= 0.002. 



Figure 33. Efficiency p(d,5) for Positive Taper. 


64 





















In Figure 34 the efficiency is plotted as a function of desynchronism d for 
negative taper rates. As the taper rate increases from 5 =0^-37t, the efficiency 
decreases from4% —> 2.6% . The peak of the efficiency curve continues to stay at small 
values of desynchronism. 



Figure 34. Efficiency Ti(d,5) for Negative Taper. 


For modest taper rates 5 =-67t ,-7;t ,-87t, the efficiency increases from 
T] = 2.9% 3.3%. With taper rate of 5 = -97t, maximum efficiencyr] =4.2% is reached 

at d = 0.002. Again sidebands are present and responsible for this enhancement in 
efficiency. For 6 =-1471 ->-2071 the peak efficiency is shifting slightly from 
d =0.002—^0.004 and efficiency diminishes from r] = 3.3%—>1.7%. Further increase 
of the taper rate 5 ^-247t, results in a decrease in efficiency to 1% at slightly larger 
d =0.006. 

Figure 35 shows the width of desynchronism Ad where the FEL can operate. 
Values of Ad against the taper rate 5 were plotted. 
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Figure 35. Effective Desynchronism Width. 


Note that A(i(5) has the same shape as the gain G(5) curve in Figure 31, with a maximum 
at 6 = 0 and decreasing rapidly as |5 1 increases. In the region where |5 1 = ±871 the rate of 
descent becomes much less rapid, indicating that something special happens there. 

The change in slope near 5 = ±87U is explained by looking at the gain curve (gain 
versus phase velocityVg) as 5 changes. Figure 36 shows that there are two gain peaks of 

equal height 5 = -87U . There are also two gain peaks when 5 = +87t, but on the opposite 

side of resonance (v = O) where (v < 0). These peaks act together, in either the 6 = ±871 

cases, to further amplify the optical field. The transition in plateau at 5 = ±871 results 
from the FEE switching for one gain peak to a second. 
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**** Gain and Phase Curves **** ) 

j=10 

fi=-87l N=41 5^=0 

0.2431 

0.0 

-0.2431 
0.16 

0.0 

-0.16 



-30 


30 


Figure 36. Comparable Peaks in Gain Spectrum, 5 =-87t. 

F. ENERGY SPREAD 

The energy spread resulting from the FEL interaction can potentially limit the 
performance of a high power FEL using energy recovery. Energy recovery is desirable 
because it can significantly increase the EEL efficiency and reduce the size of beam 
dump. In the experiment we describe, electron beam recirculation is considered feasible 
only if the induced energy spread. Ay/y is less than 6%. If it is larger, then bending 

magnets are not able to feed back the beam to the accelerator and the EEL stops 
operating. 

In order to take into account that limiting factor, we measure the width Av^ from 
the electron distribution f(y,n) in the FEL pulse evolution simulations. Accordingly we 
calculate the energy spread using Ay /y = Av^ / 4nN, where N = 41. In Eigure 37, the 
energy spread induced by the EEL interaction is plotted as a function of desynchronism d 
and positive taper 6. 
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The induced electron beam energy spread for the conventional undulator (5=0) 

has a peak energy spread of 14% occurred at d = 0.002. The strong fields cause the 
trapped-particle instability and consequently the energy spread is moved outside the 6% 
design tolerance. The induced energy spread decreases with increasing the 
desynchronism, and falls below the desired level at a desynchronism ^alue of d = 0.04, 
and remains below the 6% for all larger values of d. 

For the positive tapered undulator, the induced energy spread is slightly 
decreased. For some values of desynchronism d < 0.05 and taper 5 < IOti we are still 
outside of the limit that imposed for safe recirculation. Unfortunately this is the region in 
which we have peak efficiency. A further increase in the taper rate 5 >1071 results a 
smaller value for d for which the energy spread is acceptable. The minimum 
desynchronism that reduces the energy spread down to the desired 6% is given at the 
second column of the Table 9. 
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Figure 38 shows the final induced electron beam energy spread produced by 
inversely tapered undulators. The energy spread Ay/y, as a function of desynchronism d 
and tapering rate 5, is plotted. 



Figure 38. Electron Energy Spread for Negative Taper. 


The electron energy spread curves follow the same trend as the power curves. For 
high power, there is a large energy spread. For taper 6 = -971 the maximum value is 
18% at d =0.002. Recall that for 6 =0, the maximum value is 14% at d =0.002. 
After d = 0.05 the energy spread drops to less than 6% for all the tapering rates and 
there is no constraint to the PEL operation. The inverse tapered undulator does not appear 
to be desirable for recirculating the electron beam. This is the trade-off of the 
enhancement in the efficiency with the inverse taper. The minimum desynchronism for 
which the energy spread constraint is not met for the inverse tapered undulators is given 
at the fifth column of the Table 9. 
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Table 9 is obtained from previous Figures and helps to summarize the values of 
desynchronism d and taper 5, where we have acceptable results that meet the 
requirements. 

• min d for (Ay /y < 6%, acceptable energy spread) 

• max d for T] > 0 (FEL operation) 

In this way we are able to determine the operational range of the FEL as a function of d 
and 5. 


Negative 5 

Min d 

Max d 




-371 

0.05 

0.36 

-671 

0.05 

0.22 

-771 

0.05 

0.15 

-Stc 

0.05 

0.1 

-971 

0.05 

0.08 

-IOti 

0.05 

0.08 

-1271 

0.05 

0.08 

1471 

0.05 

0.06 

-1671 

0.03 

0.05 

-1871 

0.02 

0.04 

-2071 

0.01 

0.03 

-2271 

0.02 

0.03 

-2471 

0.02 

0.03 


Positive 5 

Min d 

Max d 

0 

0.004 

0.4 

371 

0.05 

0.36 

671 

0.05 

0.22 

771 

0.05 

0.16 

87t 

0.05 

0.1 

9k 

0.05 

0.09 

10k 

0.03 

0.09 

12k 

0.01 

0.08 

14k 

0.02 

0.06 

16k 

0.02 

0.05 

18k 

0.01 

0.04 

20k 

0.01 

0.03 

22k 

0.02 

0.03 

24k 

0.006 

0.008 


Table 9. Desynchronism Values Meeting the Energy Spread Limitation and Efficiency. 
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G. EXPERIMENTAL AND SIMULATION RESULTS -COMPARISON 

In the experiment, tapers were obtained by inserting precision shims at either end 
of the wiggler. Dial gauges on either side of the wiggler measured the position and gap of 
the wiggler for each taper value. Since the variation of the wiggler field with gap size is 
known, the field taper can be calculated from the gap taper [Ref. 40]. 

Three values of taper, 5 =-77t,0,107t were examined with the following 
experimental results compared to simulations. 


Efficiency p 



Experiment 

Simulations 

5=-77C 

2.0% 

1.54% 

5=0 

1.5% 

1.0% 

5=107t 

0.62% 

0.67% 


Energy Spread Ay/y 



Experiment 

Simulations 

5=-77C 

2.3% 

3.1% 

5=107t 

1.6% 

1.8% 


• Desynchronism Width d 



Experiment 

Simulations 

5=-ln 

0.15 

0.10 

5=0 

0.2 

0.35 

5=1071 

0.075 

0.07 


Table 10. Experiment and Simulations Results. 

We used the value of efficiency that allows safe recirculation. Eor inverse taper 6 = -Vtt , 
we disregard the highest efficiency p =3.2% at d=0.004 because it gives energy spread 
Ay /y -14% , which is almost twice the allowed limit. The same tactic was followed for 
positive taper and so we did not accept p =1.9% at d=0.002 since Ay /y -15% . Eor the 
zero taper case, the efficiency was lower than the experiment. The peak efficiency of 
p = 3.9% , at d =0.002, was accompanied by energy spread Ay/y -14% . It is possible 

to reduce the energy spread by increasing desynchronism d. hi d = 0.02 the efficiency is 
p =1.54%, while the energy spread is Ay /y = 7.5% . At d =0.03 the energy spread is 
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Ay/y = 5.7% with ri=1.32%. For all the above cases, safe recirculation near 
Ay /y ~ 6% was the goal. The FEL interaction worked well above the engineering design 
for the specific experiment. 

For the energy spread Ay /y = Av /47t we used Av from simulations, when the 
desynchronism d was at the upper limit and the efficiency was marginal. In that region 
the optical and the electron pulse fail to overlap and power was reduced to zero. For the 
inverse taper 6 =-771 that happened at d =0.1, where the efficiency was r\ =0.47% 

and for positive taper 6 = IOtt that condition was met at d = 0.05 with r] = 0.6% . 

H. CONCLUSIONS 

For small values of taper and with small values of d, the FEL exhibits the trapped- 
particle instability. Increasing either desynchronism d or taper rate 5 the EEL easily 
moves to a more stable regime, even though the efficiency is lower. This stability is 
needed in some applications, where the absence of side bands is crucial. Eor example, 
the use of the FEL in human surgery requires a narrow energy spectrum in order to 
achieve the exact desirable wavelength (A- = 6.45 |im) [Ref. 41]. 

Negative taper produces better extraction efficiency than positive taper as 
theoretical studies by Saldin, Schneidmiller and Yurkov have shown [Ref. 37]. 
Unfortunately this is not accompanied by reduction in the energy spread as p is directly 
related to Ay/y , recall p =<Aymc^ >/ymc^. At 5 =-9%, we have a 0.3% increase 
from zero taper, which results p = 4.2%. In this manner, we not only enhance the 
efficiency, but move away from the trapped-particle instability region. 

When the EEL is operated with higher Q = 50, the gain per pass is significantly 

more than the loss per pass. One result was that the range of desynchronism values when 
the EEL operates with a gain above threshold increased compared to the lower Q case. 
Eigure 39 shows the desychronism curve width Ad plotted against the positive and 
negative taper rate 6 . The superimposed G(5) shows the clear correspondence between 

increased gain G and increased desychronism curve width Ad . 
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Figure 39. Gain G(5) and Desynchronism Width Ad(5). 

The maximum gain occurs with no tapering, 5=0, as would be expected. As the 

taper rate 5 increases or decreases, the gain decreases from 160% to about 40% 
symmetrically around 5=0. At positive or negative taper rates of about 5 = ±871 , the 
gain plateaus decrease much more slowly with the increasing taper rate. While the taper 
rate increases in magnitude from |51 = Sti out to |51 = 2471 , the gain decreases from about 
40% down to about 10%. 

The gain plateau at large values of taper begins at a value of taper 5 that causes 
the gain spectrum to change shape so that there are two peaks with nearly the same peak 
gain. As the taper is increased and the peak available gain decreases, the gain spectrum 
acquires multiple peaks with comparable gain. As the primary gain peak decreases other 
surrounding gain peaks increase in comparison. As the peak available gain increases, or 
decreases, at any value of tapering, the desynchronism curve width Ad increases or 
decreases correspondingly. With more gain available there are more values of 
desynchronism that are above the threshold making the desynchronism curve wider. 

The simulations and the experimental data are in good agreement with analytic 
theory. Theory indicates that as taper increases, the desynchronism width decreases. The 
experiment shows this trend, but not as clearly as the simulation data. There is only one 
significant deviation between simulation and experimental results. The desynchronism 
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width for zero taper is A d = 0.35, compared to Ad = 0.2 in the experiment. The 

disagreement has led to reconsideration of the experimental procedure for this data point 
[Ref. 41]. 

The wiggler design can be important in the energy recovery FELs. The specific IR 
TJNAF EEL demo was designed to operate in the safe mode where Ay /y < 6% , even if 
it were capable of higher extraction efficiency Current developments at TJNAE make 
feasible the recirculation of the electron beam for all values of Ay /y that we found in 
the simulations runs. 
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VI. CONCLUSIONS 


The missile defense capabilities currently under development must not only 
review past wars but more importantly plan for the future. Tactical missile threats still in 
existence from past wars and their capabilities are rapidly improving with the result that it 
will be much more difficult to counter such threats in the future. The current development 
of weapons is centered on shortsighted decisions based on current threats and technical 
challenges. By the time some weapons are deployed they are outdated and ineffective. 
Light-based weapons offer the solution to this problem since they possess layered 
capability to counter a potential threat. 

PEL defensive capability prompted the U.S. Congress to approve a $15 million 
grant to the Thomas Jefferson National Accelerator Facility. At present, the TJNAF FEE 
is the most powerful FEE in the world and operates at more than 2 kW average power for 
many hours at a time. One transitional step for the 1 MW output power, required for a 
weapon, is the proposed 100 kW upgrade. Multimode simulations used at NFS show that 
this development is feasible. 

For the TJNAF lOOkW FEE parameters, the inverse taper undulator of A = -7U 
achieved the highest final power at a desynchronism of d = 0.04 for Q = 4.2. The output 

power produced was P~120 kW, which is above the 100 kW objective. The small 
values of the induced electron energy spread (-3.5%) made it possible to explore a design 
with 10% losses (2=10) in an optical cavity. The higher Q factor was chosen to 
increase the efficiency of the system, while maintaining the induced energy spread well 
below the 15% engineering limit. The results related with the 2 =10 design were 

• the inverse step-taper undulator A = -71 achieved the highest final power 
F = 190 kW at a desynchronism of d = 0.01. 

• The inverse linear taper 6 = -27U achieved the highest final power 
P = 200 kW at a desynchronism d = 0.01. 

The validity of our results have been compared to experiments conducted at the 
TJNAF FEE facility. Various taper rates and desynchronism values were studied. The 
purpose of the simulations was to compare experiments with theory, and, by using a 
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wider range of parameters than allowed in the experiment, to extract more physical 
meaning from the results. The simulations and the experimental data were consistent with 
analytic theory. Theory indicates that as a taper increases, the desynchronism width 
decreases. The experiment shows this trend, but not as clearly as the simulation data. 
There is only one significant discrepancy between simulation and experimental results, 
that has resulted in ongoing discussion about the validity of the experimental observation 
in this case. 

Thus it is possible to cone lude that FELs have the potential to be developed into a 
compact, reliable and efficient weapon. This development leads not only to a new 
military era but opens the door for commercial and scientific applications. The extremely 
high industrial and technological momentum that exists makes EEL development feasible 
with high level pay-offs. 
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